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Summary 


A  continuing  effort  is  being  made  to  maintain  a  closely  knit 
experimental  and  theoretical  program.  Theoretical  work  on  the  use  or 
Hermlte  functions  to  compute  Hilbert  (Kramer  Kronig)  transforms  is 
completed  and  has  been  successfully  tested  on  several  systems  and  applied 
to  the  analysis  of  optical  data  on  single  crystalline  and  amorphous  thin 
films  of  Sb2Se3. 

A  satisfactory  method  of  preparing  thin  amorphous  films  of  controlled 
composition  by  coevaporation  has  been  developed  and  improved  upon  in  the 
past  six  months.  Differential  thermal  analysis  has  yielded  Important 
information  on  the  compositional  dependence  of  the  glass  to  crystalline 
transition  in  SbxSey  amorphous  films. 

'Mossbauer  studies  have  yielded  important  information  on  the  short 
range  order  of  amorphous  films  and  on  the  structural  chemistry  of  mixed 
chalcogenide  systems.  This  technique,  where  applicable,  appears  to  fulfill 
a  need  for  analysis  on  a  microscale  where  more  conventional  analytical 
tools  such  as  X-ray  and  electron  microscopy  are  inadequate.  However, 
initial  studies  of  SbxSey  amorphous  films  by  electron  microscopy  are 
yielding  very  interesting  data  but  a  full  interpretation  of  the  data 
awaits  compositional  analysis  of  the  films. 

-In  the  optical  studies  of  thin  amorphous  films,  one  of  the  mere 
significant  findings  has  been  the  close  correspondence  between  the  energy 
band  gaps  in  the  crystalline  and  amorphous  form  of  Sb2Se3  and  the  apparent 
lack  of  band  tailing.  Further  optical  measurements  on  single  crystal  Sb2Se3 
in  the  neighborhood  of  the  band  edge  as  a  function  of  crystal  orientation 
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have  supported  the  initial  findings  of  a  direct  transition.  We  are  currently 
pressing  ahead  with  resistivity  versus  temperature  measurements  on  both  the 
amorphous  and  single  crystal  forms  to  obtain  thermal  activation  energies 
for  comparison  with  the  optical  values. 

Several  different  approaches  to  modulated  reflectivity  (and  modulated 
absorption)  are  under  development.  This  data  should  prove  invaluable  in 
the  analysis  of  the  electronic  structure  of  the  antimony  chalcogenides. 

Along  the  same  lines,  the  setting-up  of  appartus  for  the  very  powerful 
technique  of  photoemission  is  at  a  final  stage  of  completion  and  is  already 
yielding  results  on  the  density  of  states. 

As  evidenced  by  the  large  number  of  references  to  Bull.  Am.  Phys. 

Soc.  much  of  the  work  is  still  In  the  initial  stages.  In  keeping  with 
past  practice  we  have  Included  preprints  of  work  currently  being  submitted 
for  publication  to  elaborate,  for  the  interested  reader,  on  details  not 
suitable  for  this  general  progress  report.  A  summary  of  the  more  signi¬ 
ficant  findings  of  this  program  to  date  will  be  presented  at  the  forthcoming 
International  Symposium  on  Amorphous  and  Liquid  Semiconductors  and  an 
abstract  is  appended  to  this  report.  Reprints  of  papers,  published  since 
the  last  six-monthly  report,  are  also  appended. 


THEORETICAL  RESEARCH 


I.  The  Kramers-Kroniq  (Hilbert)  Transform 

During  the  last  six  months  the  work  initiated  earlier  on  the  use 
of  Hermite  functions  to  compute  Hilbert  (Kramers-Kronig)  Transforms  has 
been  completed.  A  paper  outlining  the  general  mathematical  technique 
has  been  submitted  for  publication  and  a  preprint  was  included  in  the 
previous  technical  report.  This  paper  concentrates  upon  the  application 
of  the  method  to  functions  of  compact  support  such  as  densities  of  states. 
Hilbert  Transforms  of  such  functions  are  useful,  e.g.,  in  the  computation 
of  optical  functions  from  the  results  of  electronic  structure  calcula¬ 
tions.  Briefly,  a  convergent,  matrix  representation  of  the  Hilbert 
Transform  kernel  has  been  employed  to  compute  Hilbert  Transforms  of  two 
densities  of  states,  1)  the  classic  simple  cubic  s-band  (triple  cosine) 
density  of  states  and  2)  the  electronic  density  of  states  of  fee  palla¬ 
dium1.  The  results  indicate  that  critical  point  structure  does  not  shift 
In  frequency  under  a  proper  Hilbert  Transform.  Details  of  the  computat 
tional  method  and  its  analytical  development  are  found  in  the  paper. 

The  application  of  the  above  technique  to  optical  data  analysis,  i.e. 
to  the  extraction  of  optical  functions  from  experimental  near-normal  in¬ 
cidence  reflectivity  measurements  has  also  come  to  fruition.  A  preprint 
of  a  paper  which  is  now  in  its  final  stages  of  preparation  is  attached 
to  this  report  as  Appendix  A.  A  preliminary  report  of  the  results  of  the 
technique  as  applied  to  Sb2Se^  crystalline  samples  was  presented  at  the 
March  Meeting  of  the  American  Physical  Society  in  Cleveland  ^(Appendix  B). 
Further  application  of  the  technique  has  been  made  to  thin  amorphous  films 
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of  the  composition  Sb2Se3  and  this  work  was  reported  at  the  April  Meeting 
of  the  American  Physical  Society  in  Washington^Appendix  C).  There  are 
indications  of  substantial  differences  in  the  structural  details  of  the 
optical  functions  of  crystalline  and  amorphous  Sb2Se3-  Further  investi¬ 
gations  in  this  area  are  underway.  Basically  the  method  employs  the  re¬ 
flectivity  of  systems  of  simple  charged  oscillators  to  extrapolate  the 
experimental  data  outside  the  measured  data  range  and  to  convert  the 
function  whose  Hilbert  Transform  needs  to  be  computed  to  one  of  compact 
support.  The  methods  employed  in  the  previous  work  described  above  are 

then  applicable.  As  a  check  upon  the  power  of  the  technique  it  has  been 
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applied  to  reflectivity  data  on  BaF2  in  the  far  ultraviolet.  This  data 
has  a  wealth  of  sharp  structure  arising,  most  likely,  from  deep  excitons 
and  provides  a  test  of  the  resolution  of  the  computational  technique. 

II.  Electronic  Structure 

Some  progress  has  been  made  on  the  computation  of  the  electronic 
spectra  of  amorphous  materials.  The  method  we  have  explored  is  similar 

c 

in  spirit  to  the  work  of  Maschke  and  Thomas  for  electronic  structure 
in  amorphous  solids  with  short-range  order.  We  are  planning  tc  apply 
our  method  to  specific  amorphous  solids.  A  computer  program  to  compute 
electronic  energies  of  a  model  of  amorphous  Ge  and  the  associated  density 
of  states,  is  under  development.  The  computational  problems  are  severe 
both  as  to  complexity  and  as  to  computational  time  required.  The  problem 
is,  however,  tractable  and  with  continuing  work  will  be  completed. 

In  another  vein  work  has  been  initiated  to  attempt  to  make  use  of 
the  results  recently  obtained  by  Romeris®  to  connect  the  electronic 
structure  of  amorphous  so^ds  to  that  of, crystals  of  similar  composition. 
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Romeris  has  shown  that  critical  point  structure  behaves  in  a  particular 
manner  as  a  "crystal"  potential  varies  from  being  truly  periodic  to  being 
almost  periodic.  Since  there  is  evidence^  of  critical  point  structure  in 
the  optical  functions  found  for  amorphous  films  of  Sb2Se3  it  may  be  feasi¬ 
ble  to  correlate  this  structure  in  a  definite  manner  with  that  found  in 
crystalline  Sb2Se^.  We  have,  therefore  begun  to  explore  methods  for 
tracking  the  behavior  of  critical  point  structure  with  disorder. 


4 


MATERIALS  RESEARCH 


I.  Amorphous  Thin  Films 

A  large  number  of  amorphous  films  of  various  compositions  in  the 
Sb-Se  system  have  now  been  prepared  by  coevaporation.  The  compositions 
were  determined  by  micro-probe  analysis  in  conjunction  with  the  Geo¬ 
physical  Laboratory  of  the  University  of  Chicago.  A  fair  degree  of 
reproducibility  has  been  obtained  but,  because  of  the  critical  depen¬ 
dence  of  optical  properties  on  composition7  (Appendix  B)  the  deposition 
controllers  (Sloan  Omni  IIA)  have  had  to  be  upgraded.  The  most  signi¬ 
ficant  improvements  are:  newly-developed  circuit  boards  have  been  in¬ 
stalled,  gold-plated  quartz  crystals  of  circular  geometry  are  now  employed 
and,  most  importantly,  the  frequency  and  deposition  rate  of  each  source 
is  continuously  monitored  externally  by  chart  recorders  throughout  depo¬ 
sition.  Despite  these  improvements  control  of  stochiometry  to  within  1 
atomic  percent  cannot  be  guaranteed. 

In  an  attempt  to  obtain  even  closer  control  of  composition  two  further 
approaches  will  be  tried:  1)  sputtering,  under  equilibrium  conditions: 

2)  coevaporation  of  a  graded  composition  across  an  extended  substrate. 

Both  systems  are  under  construction,  the  former  being  nearly  completed. 

II.  Crystalline  Thin  Films 

The  three- temperature  technique®  was  successfully  applied  to  the 
deposition  of  a  stochiometric,  polycrystalline  film  of  Sb2Se3«  There 
were  difficulties  in  achieving  adequate  deposition  rates  which  are 
determined  by  both  substrate  temperature  and  melt  temperature.  Too  fast 
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a  deposition  rate  produced  a  film  with  small,  irregular  crystal  grains. 
With  too  slow  a  rate  there  is  virtually  no  deposition.  It  was  found 
that  a  fused-quartz  substrate  temperature  of  450°C  and  a  melt  temperature 
of  630°C  produced  a  satisfactory  film  for  optical  studies. 

III.  Single  Crystals 

A  new  single  crystal  growth  apparatus  is  under  development  using  a 
tin-oxide  coated  quartz  tube  for  both  the  background  heating  elements  and 
the  zone  heating  element.  These  tubes  have  the  advantage  that  they  will 
develop  temperatures  well  above  600°C  while  remaining  very  transparent 
to  allow  good  visibility  of  the  molten  zone  and  hence  good  control  of  the 
growth  process. 

IV.  Differential  Thermal  Analysis 

DTA  measurements  have  been  carried  out  on  various  compositions  of 
Sb-Se  films  in  the  temperature  range  of  30°C  to  300°C9  (Appendix  C).  The 
glass  to  crystalline  phase  transitions  of  the  samples  were  investigated. 
This  transition  was  found  to  be  endothermic  and  located  ar  ^190°C  for 
Sbj)Se3.  Those  films  that  were  on  the  Se-rich  side  of  SbgSeg  showed 
amorphous-crystalline  transition  for  Se,  which  was  located  at  ,vl20oC,  as 
well  as  the  Sb2Se3  transition.  Plans  have  been  made  to  incorporate  a 
more  sensitive  arrangement  of  thermo-couples  more  suited  to  a  thin  film 
geometry. 


STRUCTURAL  INVESTIGATIONS 


I.  Mossbauer  Studies 

Work  performed  in  collaboration  with  Dr.  S.  Ruby  of  Argonne  National 
Laboratory  produced  the  very  significant  finding  that  all  Sb  atoms  in 
SbxSe^  (Including  Sb2Se3 )  amorphous  solids  sit  in  closely  similar  environ¬ 
ments  in  contrast  to  the  Sb2Se.j  crystal  where  there  are  two  distinct 
sites10  'Appendix  B).  This  is  of  particular  significance  in  that  it  lends 
support  to  the  view  that  the  evaporated  fl ims  are  in  reality  amorphous 
and  not  "micro-crystalline",  since  we  now  know  that  the  short-range- 
order  differs  in  the  two  forms.  This  work  will  be  submitted  for  publi¬ 
cation  and  a  preprint  is  attached  (Appendix  D). 

In  order  to  investigate  cross-linking  effects  of  Group  IV  elements 

in  amorphous  Group  V  chalcogenides,  Mossbauer  studies  using  ^Sn  in 

a  BaSn03  source,  have  been  carried  out  on  Snx(As2Se2)lx  glasses  with 

varying  concentrations  of  tin,  i.e.  (x  =  2.5,  5.0,  8.0,  10.0  At.  %)  and 

on  the  two  crystalline  samples  of  SnSe2  at  liquid  nitrogen  and  liquid  He 

temperatures! 1  (Appendix  B).  Sn  (AsoSeo),  glasses  samples  were  prepared 

x  c  0  '“X  0 

by  quenching  the  melts  into  a  brine  bath  at  200K. 

The  limit  of  tin  solubility  in  the  quenched  glassy  state  is  tenta¬ 
tively  set  by  Kolomiets  et  al1^  at  7  At.  %.  Thus  8.0  and  10.0  At.  %  Sn, 
glass  samples  may  have  a  crystalline  component.  X-ray  patterns  for  all 
the  samples,  however,  were  similar  and  characteristic  of  amorphous  systems. 
The  major  phase  present  in  a  deliberately  crystallized  AsgSes  +  10  At.  %  Sn 
sample  is  SnSe  and  the  corresponding  amorphous  sample  also  contains  a 
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small  amount  of  SnSe.  A  critical  Mossbauer  spectra  for  10  At.  %  Sn 
amorphous  sample  at  4.2°K  shows  an  increase  of  absorption  by  about  15% 
compared  to  that  at  liquid  nitrogen  temperature.  The  overall  character 

1  O 

of  the  spectra  is  similar  to  that  obtained  ty  Borisova  et.  al.  ,  the 
major  difference  occurred  in  the  intensities  of  the  spectra  correspond¬ 
ing  to  SnSe2  and  SnSe.  The  average  shift  in  all  systems  corresponding 
to  SnSe2  is  1.68  mm/sec  which  is  slightly  less  than  that  of  a-Sn.  On 
the  basis  of  Ruby's^  atomic  calculations,  the  electronic  configuration 
is  0.7s,  3.3p.  The  average  shift  corresponding  to  SnSe  is  3.45  mm/sec 
and  the  electronic  configuration  turns  out  as  1.6s,  2.4p. 

Thus,  it  is  clear  from  the  present  Mossbauer  investigations  of  the 
amorphous  systems  that  the  degree  of  local  order  is  very  high,  i.e., 
molecules  or  microcrystals  of  small  enough  dimensions  to  escape  x-ray 

detection  are  held  in  a  disordered  matrix.  The  microchemistry  does 

15 

appear  to  correspond  to  the  suggested  chemistry  of  Shkolnikov  .  The 
Mossbauer  spectra  for  Sn  in  the  glasses  corresponds  very  closely  to 
those  of  crystalline  SnSe  and  SnSe2*  However  the  range  of  chemical  order 
is  too  small  to  be  detected  by  x-ray  diffraction. 

II.  Electron  Microscopy 

Electron  microscopy  and  electron  diffraction  studies  were  made  on 

thin  (about  250°  8  thick)  films  of  Sb  Se  „  Amorphous  films  appear  granu- 

*  y 
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lar  with  irregularly  shaped  particles  of  about  1  x  10  cm  diameter.  The 
corresponding  diffraction  pattern  consists  of  three  diffuse  bands.  One 
can  observe  crystal  growth  from  these  amorphous  films  as  caused  by  the 
heat  of  the  electron  beam.  Depending  on  the  sample  composition,  diff¬ 
erent  types  of  ordering  are  observed:  1)  bizarre  looking  dondri te-1 ike 
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needles  of  about  1000  a  width  with  a  complicated  fine  structure,  2)  fibres 
of  about  60  8  width  which  are  lined  up  parallel  to  one  another  or  are 
spiralling  around  each  other,  3)  rounded  domains  of  micron  size  with  narrow 
subdomains  of  about  100  x  1000  8  which  have  different  alignments  for 
different  domains. 

The  diffraction  patterns  of  sections  50  u  in  diameter  of  the 
crystallized  films  range  from  spotty  powder  patterns  and  fibre  patterns 
to  single  crystal  patterns.  Typical  is  a  one  dimensional  ordering 
resulting  in  irregularly  spaced  spots  on  equally  spaced  lines.  The 
ordered  Jimension  varies  somewhat  around  3.24  8.  It  is  not  known  at 
this  date  if  the  corresponding  reflections  relate  to  the  "forbidden" 

(100)  and  (010)  reflections  from  Sb203  structure  or  belong  to  a  completely 
new  crystal  type.  Definitely  not  compatible  with  the  Sb203  structure  are 
recurring  hexagonal  patterns  which  cannot  be  attributed  either  to  selemium 
or  antimony. 

A  complete  interpretation  of  the  results  will  be  attempted  when 
all  the  sample  compositions  have  been  analyzed. 


JL 


2 


9 


OPTICAL  PROPERTIES 


I.  Transmission  and  Reflectivity  Measurements 

Several  different  methods  have  been  employed  to  determine  the 
optical  constants  n  and  k,  from  the  reflectivity  (R)  and  transmission  (T) 
of  a  thin  film  in  the  spectral  region  near  and  below  the  absorption  edge. 
The  first  technique  tried  was  to  search  the  n-k  plane  for  points  that 
would  fit  the  observed  R  and  T.  This  method  had  limited  success  because 
for  a  given  R  and  T  there  are  multiple  solutions  in  the  n-k  plane  and  the 
computer  program  was  taking  excessive  time.  The  other  method  tried  was 
that  of  quadratic  interpolation.  This  metnod  is  not  yet  fully  completed, 
but  it  gives  greater  accuracy  in  less  computer  time. 

Absorption  constant  (a)  vs.  photon  energy  graphs  were  determined  by 
transmission  and  reflectivity  data  for  a  wide  range  of  compositions 
of  Sb-Se  amorphous  films7  (Appendix  B).  A  close  correspondence  was 
noted  between  the  position  of  the  optical  gap  of  the  amorphous  films 
and  single  crystal  Sb2Seg.  The  a  values  were  determined  from  106cm_1 
to  10  cm-1  for  selected  samples  using  a  wide  variation  of  sample  thick¬ 
nesses.  Thin  films  (^5  )  of  Sb0Se,  must  still  be  made  and  measured  inn 

m  Z  3 

order  to  determine  the  type  of  transition  that  occurs  (direct  o t  phonon 
assisted). 

The  reflectivity  and  transmission  of  single  crystals  of  Sb2Se3 
grown  in  a  horizontal  zone  refiner  were  measured  in  the  range  0.03  eV 
to  11  eV1®  (Appendix  B).  R  and  T  were  measured  as  a  function  of  polari¬ 
zation  in  the  cleavage  plane.  The  reflectivity  structure  was  shown  to 
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agree  with  Soblov  et.  al.17  for  the  two  polarization  directions  in  the 
cleavage  plane  in  the  energy  range  1  eV  to  4  eV;  and  for  unpolarized 
light  in  the  range  4-11  eV.  The  energy  range  0.03  to  0.5  eV  was 
covered  on  the  Beckman  IR-12  spectrophotometer;  0.5-4  eV  on  a  Cary  14  R 
spectrophotometer;  and  4-11  eV  on  a  McPherson  vacuum  U.V.  spectrophoto¬ 
meter  at  Argonne  National  Laboratory  (kindly  loaned  by  Dr.  P.  Yuster). 

The  absorption  coefficient  was  calculated  from  this  data  near  the  ab¬ 
sorption  edge  at  about  1  eV.  Different  powers  of  the  absorption  co¬ 
efficient  were  plotted  versus  wavelength  to  obtain  the  type  of  transition 
responsible  for  the  absorption  edge.  A  direct  allowed  transition  was 
found  for  both  directions  of  polarization. 

Initial  low  light  level  measurements  have  been  made  using  an  S.S.R. 
photon  counter.  Great  lengths  have  been  gone  to  in  order  to  prevent 
light  leaks  in  the  system.  The  phototubes  can  be  cooled  down  to  solid 
CO2  temperatures  for  operation  in  a  low  noise  mode.  This  arrangement 
should  prove  to  be  quite  valuable  for  the  measurement  of  optical  trans¬ 
mission  well  above  the  absorption  edge  and  also  for  small  signal  modu¬ 
lation  experiments. 

II .  Modulated  Reflectivity 

Thermally  modulated  reflectivity  measurements  were  performed  on 
single  crystals  of  Sb2Te3  with  the  light  incident  on  the  cleavage  plane 
(El Cj^Appendix  C).  The  data  was  taken  in  the  visible  and  near  infra¬ 
red  regions  and  was  found  to  agree  with  critical  points  observed  in  the 
transmission  spectrum  of  thin  cleaved  sample  using  a  photon  counter. 

Electro-reflectance  apparatus  is  now  completed  which  will  be  used 
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to  measure  the  electro-modulated  reflectivity  spectrum  of  Sb2Se3  single 

crystals,  in  all  three  orientations  of  the  electric  vector,  and  amorphous 

Sb  Se  films. 

*  y 

Also  under  construction  is  apparatus  for  wavelength  modulated  re¬ 
flectivity  which  involves  halving  the  optical  beam  from  the  exit  slit 
of  a  monochrometer  and  comparing  the  reflectivity  of  the  sample  to  light 
from  the  two  halves  by  means  of  a  lock-in  amplifier. 

Construction  of  apparatus  for  rotoreflectance  is  also  underway. 
This  technique  utilizes  a  rotating  polarizer  to  measure  the  changes  in 
reflectivity  as  a  function  of  electric  vector  orientation. 

III.  Photoconductivity 

In  the  photoconductivity  area  most  work  was  concentrated  on  study 
of  amorphous  SbxSey.  The  spectral  dependence  of  photoresponse  and 
photoconductive  gain  depend  critically  on  sample  composition.  The 
photoconductivity  activation  energy  (as  determined  by  Moss'19  criterion) 
generally  agrees  with  the  optical  activation  energy  as  determined  from 
transmission  measurements^  (Appendix  B)  in  case  of  stochiometric  amor¬ 
phous  Sb2Se3  being  equal  to  1.1  eV,  which  is  very  close  to  that  of  single 
crystalline  Sb2Se3-  We  feel  that  this  point  is  of  fundamental  importance 
as  it  shows  that  the  gap  is  approximately  the  same  in  amorphous  and 
crystalline  Sb^Sej;  contrary  to  previously  published  work  by  Kolomiets 
et.  al.^0  whose  results  were  probably  affected  by  deviations  from  stochio 
metry.  The  effect  of  trapping  states  inside  the  forbidden  gap  on  photo- 
conductive  gain  and  on  the  response  time  was  investigated  using  a  dc 
light  "bias"  of  variable  intensity  to  fill  traps  to  any  desired  degree. 

In  order  to  provide  meaningful  comparison  of  the  photoconductive 
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properties  of  amorphous  Sb2Seg  with  single  crystalline  Sb2Se2»  a 
detailed  investigation  of  the  photoconductive  properties  of  the  single 
crystalline  phase  was  performed?'  (Appendix  B).  The  onset  of  photo¬ 
conductivity  was  found  to  depend  on  orientation  of  electrical  vector  E 
of  polarized  light  with  respect  to  three  principal  crystallographic 
axes  in  agreement  with  dependence  of  optical  absorption  on  the  above 
parameters. 

An  interesting  feature  of  the  photoconductive  kinetics  of  SbgSe^  is 
a  rather  long  response  time  -  even  in  the  region  of  fundamental  absorption 
the  response  time  is  of  the  order  of  250-300  micro-seconds,  which  can  be 
explained  only  by  strong  trapping  effects  which  obscure  the  surface  re¬ 
combination.  This  long  response  time  probably  also  gives  rise  to  photo¬ 
sensitivity  over  wide  spectral  range  (1-3  eV). 

So  far,  we  are  able  to  make  the  following  conclusions:  The  spectral 
response,  particularly  the  activation  energy  of  amorphous  stochiometric 
Sb2Se3  is  very  similar  to  that  of  single  crystal  Sb2Se.j.  The  onset  of 
photoconductivity  in  amorphous  Sb2Se2  is  at  least  that  sharp  as  in  the 
case  of  the  single  crystal  modification,  a  fact  which  seems  to  be  contra- 
dictory  to  effects  of  band  tailing  as  found  by  Weiser  and  Brodsky4,4,  for 
amorphous  As2Te3.  Also  the  response  time  and  photosensitivity  is  at 
least  two  orders  of  magnitude  lower  than  that  in  the  crystalline  phase 
indicating,  together  with  the  previous  findings,  that  trapping  states  at 
band  edges  in  amorphous  phase  of  SbgSe^  seem  to  be  electrically  much  less 
effective  than  in  the  single  crystal  phase. 

IV.  Photoemission 

Construction  of  a  ultrahigh  vacuum  (10"^  -  10“^  torr)  phototube 
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has  been  completed  together  with  provision  for  sample  cleaving  in  situ. 
Presently,  an  attachment  is  under  construction  for  the  evaporation  of 
amorphous  Sb-Se  films  in  situ. 

Also  completed  is  the  whole  detection  and  electronic  circuit  em¬ 
ploying  a  retarding  field  detector.  A  few  preliminary  runs  show  that 
the  apparatus,  i.e.,  McPherson  225  Vacuum  UV  Monochromator,  all  vacuum 
parts  and  all  electronic  equipment,  is  in  good  working  condition  and 
initial  density  of  states  curves  have  been  recorded.  These  curves  will 
be  compared  with  density  of  states  from  band  structure  calculations  for 
amorphous  and  crystalline  Sb2Se3.  In  order  to  obtain  better  resolution 
of  transitions  at  critical  points  in  Brillouen  zone  partially  obscured 
by  phonon  broadening  and  by  the  absoptive  background,  first  and  second 
derivatives  of  the  density  of  states  vs.  energy  will  be  also  recorded 
simultaneously.  For  this  purpose  doublers  and  triplers  of  the  modula¬ 
tion  frequency  have  been  constructed. 
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TRANSPORT  PROPERTIES 


I.  Resistivity  vs.  Temperature 

Work  has  progressed  In  the  low  temperature  measurement  of  very 
high  resistance  (^10^  ohms)  thin  films.  These  measurements  have  been 
carried  out  In  a  "Cryotip",  which  uses  the  Joul e-Thompson  expansion  of 
He  and  H2  to  provide  temperatures  from  ambient  to  4.2°K.  A  four-probe 
configuration  of  evaporated  gold  contacts  Is  used  to  minimize  contact 
resistance.  The  main  problems  encountered  thus  far  are  those  of  vacuum 
Integrity,  shielding  of  circuits  and  ground  loops  between  circuits. 
These  have  been  solved  to  the  extent  that  a  standard  resistance  of  the 
same  order. of  magnitude  as  the  film  has  been  measured  In  the  apparatus. 
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ABSTRACT 

Herein  a  new  technique  for  performing  Kramers-Kronig  (Hilbert) 
transforms  upon  experimental  normal  incidence  reflectivity  data 
is  presented.  The  method  is  based  upon  development  of  the 
Hilbert  principal  value  kernel  and  suitably  defined  functions 
closely  related  to  the  data  in  the  space  of  Hermite  functions. 

Well  defined  procedures  for  extrapolating  the  data  into  unmeasured 
regions  of  the  spectrum  are  employed. 

The  method  is  applied  to  reflectivity  data  taken  on  two 
solids  and  the  results  are  analyzed  as  to  their  reliability  as 
is  the  general  utility  of  the  technique.  An  important  conclusion 
drawn  from  the  results  is  the  invariance  of  the  location  of 
critical  point  structure  under  the  Kramers-Kronig  transform. 

Some  of  the  details  of  the  numerical  procedures  and  their  con¬ 
vergence  are  also  briefly  discussed. 


APPENDIX  A 

-  3  - 


I.  INTRODUCTION 


Great  progress  has  been  made  in  extending  the  frequency 
region  over  which  the  optical  spectroscopy  of  solids  is  regularly 
performed.  It  is  presently  relatively  commonplace  for  measure¬ 
ments  of,  e.g.,  near-normal  incidence  reflectivity,  to  be  carried 
out  from  far  infra-red  through  vacuum  ultra-violet  frequencies. 
Reflectivity  measurements  on  bulk  single  crystals  over  wide 
spectral  ranges  are  made  especially  valuable  by  the  possibility 
of  performing  a  Kramers-Kronig  analysis  upon  the  data  to  obtain 
the  phase  of  the  complex  reflectivity  and  hence  all  of  the  optical 

functions  of  a  given  material.  This  general  procedure  has  been 

/  j-'l ) 

adopted  by  many  workers'  •  and  several  computational  techniques 
have  been  developed  for  carrying  out  the  requisite  integral 
transforms.  In  Section  II  these  are  briefly  and  critically 
reviewed.  Also  in  that  section  the  method  we  have  employed  in 
this  work  is  outlined.  The  general  procedure  has  been  separately 
detailed  in  a  somewhat  different  context'  '  and  so  the  intent 
here  is  to  concentrate  upon  those  aspects  of  Kramers-Kronig 
analysis  and  of  our  method  which  are  specific  to  the  analysis 
of  optical  data. 

In  Section  III  the  results  of  the  application  of  this  new 

(  \ 

method  to  normal  incidence  reflectivity  data  on  SbgSe^  '  and 

(  7  ) 

BaPg'  '  single  crystals  are  presented.  The  reflectivity  data 
In  these  two  cases  were  measured  over  different  spectral  regions 
and,  additionally,  the  two  sets  of  data  are  quite  different  in 
structural  detail.  They  thus  provide  tests  of  somewhat  different 
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character  (e.g.,  as  to  resolution)  for  the  technique. 

In  Section  IV  an  analysis  of  the  accuracy,  reliability, 
and  systematic  use  of  the  method  is  given.  Studies  of  the 
convergence  of  the  necessary  numerical  expansion  procedures  are 
presented  as^a  discussion  of  some  general  conclusions  on  the 
characteristics  of  the  Hilbert  transform  as  it  relates  to  the 
optical  properties  of  solids.  Possible  applications  to  other 
problems  in  solid  state  physics  are  considered. 


II.  THE  KRAMERS-KRONIO  (HILBERT)  TRANSFORM 
APPLIED  TO  OPTICAL  DATA 


If  it  were  possible  to  measure  the  (near- normal  incidence) 
fractional  coefficient  of  reflected  Intensity,  or  reflectance, 

R(  u>) »  of  a  bulk  solid  sample  (for  a  given  orientation  and  polari¬ 
zation)  in  vacuo  (air)  for  all  frequencies  then  the  phase  of  the 
complex  (amplitude)  reflection  coefficient  r(  <„)  =  o))e1(p(<«) 

(  8  ) 

could  be  determined  from  the  Kramers-Kronig  dispersion  relation'  ' 


l/2 

where  r(u)l)  =  R  '  ( u,*  )*,  and  R(  m)  has  been  extended  to  negative 

*Most  often  dispersion  relations  are  derived  for  the  real  and 
Imaginary  parts  of  comojp^  qanpaji  pppponse  functlons/'whlcfc  of  a 
real  variable ^must  be  of  class  L‘c(-eo,  ®) .  However  such  relations 
can  also  be  derived  for  the  phase  and  logarithm  of  the  modulus 
of  complex  functions  for  which  £lm  r(  uj)  -♦  0  at  least  as  rapidly  as 

U)-*» 

co"  where  n  >  1.  Thus  although  the  logarithm  of  the  reflectivity 
is  not  L2(-®,  +»)  a  dispersion  relation  exists  between  it  and  the 
phase.  For  causal  response  functions  the  real  and  imaginary  parts 
are  each  of  class  L2(-®,  »).  An  excellent  discussion  of  these 
points  is  given  in  Ref.  (  c  ) . 
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frequencies  as  an  even  function.  Knowledge  of  R(u>)  and  is 

tantamount  to  knowledge  of  all  the  optical  functions  since 
6(  a,)  =  n(t0)  +  ix(  cu)  is  related  to  the  complex  reflection  coefficient 
by 

fi(»)  -  (1  +  £(»))/(!  -  £(»))  .  (2) 

a  simple  homographic  transformation.  Given  r\(u>);  e(u>)>  the  com¬ 
plex  dielectric  function,  and  the  energy  loss  function, 

]»  are  easily  obtained. 

There  are  two  paramount  difficulties  in  Eq.  (1)  which  must 
be  overcome  in  order  that  this  data  reduction  be  effected.  In 
the  first  place  data  is  never  available  for  all  frequencies  due 
to  the  finite  frequency  cutoff  of  real  spectrometers.  Secondly 
the  kernel  of  the  Hilbert  transform  in  Eq.  (1)  is  highly  singular 
and,  in  fact,  that  equation  is  most  properly  regarded  as  a  functional 
and  the  kernel  as  a  generalized  function.  A  casual  numerical 
attack  upon  it  or  any  of  its  equivalents  can  be  rather  precarious 
particularly  in  its  effect  upon  the  critical  point  structure 
exhibited  by  the  spectrum  under  consideration. 

It  Is  useful  to  review  some  of  the  methods  which  have  been 
employed  to  deal  with  the  aforementioned  problems.  Generally 
the  first  problem,  that  of  finite  data  ranges  has  been  handled 
by  extrapolating  the  data  Into  unmeasured  regions.  Above  the 
highest  frequency  the  data  is  usually  extended  using  one  or 
more  terms  of  a  Lavrent  series.  Unless  great  care  is  taken 
this  may  introduce  spurious  structure  into  the  spectrum.  Below 
the  lowest  frequency  measured  the  data  is  normally  extrapolated 
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be  overcome  in  order  that  this  data  reduction  be  effected.  In 

the  first  place  data  is  never  available  for  all  frequencies  due 

to  the  finite  frequency  cutoff  of  real  spectrometers.  Secondly 

the  kernel  of  the  Hilbert  transform  In  Eq.  (1)  is  highly  singular 

and,  in  fact,  that  equation  is  most  properly  regarded  as  a  functional 

and  the  kernel  as  a  generalized  function.  A  casual  numerical 

attack  upon  it  or  any  of  its  equivalents  can  be  rather  precarious 

particularly  In  its  effect  upon  the  critical  point  structure 

exhibited  by  the  spectrum  under  consideration. 

It  is  useful  to  review  some  of  the  methods  which  have  been 

employed  to  deal  with  the  aforementioned  problems.  Generally 

the  first  problem,  that  of  finite  data  ranges  has  been  handled 

by  extrapolating  the  data  into  unmeasured  regions.  Above  the 

highest  frequency  the  data  is  usually  extended  using  one  or 

more  terms  of  a  Lavrent  series.  Unless  great  care  is  taken 

this  may  introduce  spurious  structure  into  the  spectrum.  Below 

the  lowest  frequency  measured  the  data  is  normally  extrapolated 
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using  a  Taylor  series.  Knowledge  of  the  static  dielectric  con¬ 
stant  of  the  material  being  studied  can  make  this  a  most  reasonable 
procedure.  The  singular  kernel  in  (1)  is  usually  handled  by 
transforming  (1)  into  the  equally  valid  form 

«p( a)  Ss  f  ->*|  -  IgaLal  a,.  .  (3) 

The  inclusion  of  a  small  additional  part  in  the  denominator  of 

p 

(3),  Tl  ,  is  then  justified  by  the  principal  value  character  of 
(1)  and  if  it  is  chosen  sufficiently  small  the  finite  arithmetic 
character  of  digital  computation  allows  the  transfoim  to  carried 
out  in  an  approximate  manner.  It  is  interesting  to  examine  in 
more  detail  the  methods  used  by  some  of  the  workers  who  have 
attacked  the  Kramers-Kronig  transform  ir.  optical  problems^  tws r 
In  another  publication^  *  ^  (hereafter  denoted  I)  the 

I  *  t* 

suitably  of  the  Hermite  functions  (harmonic  oscillator  eigenfunc¬ 
tions)  as  a  basis  set  for  accomplishing  Hilbert  transforms  has 
been  demonstrated.  The  interest  there  was  in  the  Hilbert  transforms 
of  electron  (or  phonon)  densities  of  states  and  thus  the  functions 
to  be  transformed  were  of  compact  support,  i.e.,  they  vanished 
outside  some  finite  range  on  the  real  line,  and  were  bounded. 

In  the  application  of  Kramers-Krunig  analysis  to  reflectivity 
data  these  characteristics  are  not  present  however  as  will  be 
shown  suitable  treatment  of  the  data  brings  the  problem  into  a 
form  amenable  to  the  methods  described  in  I. 

In  I  it  is  shown  that  if  a  function  <x(  co)  have  a  convergent 


Hermite  expansion 
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°(“)  s  d  “n  ♦n^  (4) 

where  the  ♦,,(«■»)  are  the  Hermite  functions,  then  Its  Hilbert  transform 
partner  8(«)  also  has  a  Hermite  expansion 

eU)  a  S0n  *n(«)  (5) 

n 

where 

®n  =  m  ^nm  “m  » 

being  a  compact  matrix  representation  of  the  Hilbert  transform 

* 

"operator"  in  the  space  of  Hermite  functions.*  vanishes  unless 
n  +  m  is  odd  and,  for  integer  h,  may  be  written 

snm  =  (’Dh  +  1  h  =(m-n-l)/2  (7) 

Computation  of  the  matrix  elements  S^,  is  most  easily  accomplished 
as  in  I;  from 

+“ 

Snm  =  _mJ  dt  SGN(t)*n(t)*m(t)  .  (8) 


*  It  can  be  shown  (8)  that  causal  complex  response  functions  which  are 
the  simplest  class  of  functions  whose  real  and  imaginary  parts  are 
Hilbert  transform  pairs  are  necessarily  square  integrable  as  are  their 
real  and  imaginary  parts  separately. (8)  However  the  class  of  square 
integrable  functions  does  not  exhaust  all  possible  Hilbert  transform* 
pairs  (consider  <j>(w)  and  tnr(u)  in  (1)).  The  Hermite  decomposition 
of  the  Hilbert  transform  can  be  directly  applied  only  to  the  square- 
integrable  case  because  of  the  need  for  convergent  Hermite  expansions 
to  exist. 


30 


APPENDIX  A 
-  6  - 


The  matrix  elements  are  then  computed  by  using  the  Hermite  inte¬ 
gration  technique,  as  detailed  in  (i). 

Direct  application  of  this  technique  to  (l)  is  obviated 
since  as  (u  •*  ±  o  r(  uj)  -♦  0,  hence  ^nr(  u,)  -*  -®  as  m  -♦  ±  ®. 

The  function  is  therefore  not  a  square  integrable  function 

and  has  no  convergent  expansion  in  Hermite  functions.  A  simple 
procedure  however  removes  this  difficulty  allowing  the  full 
power  of  the  Hermite  analysis  to  bear  upon  the  structure  in  the 
optical  problem  as  well  as  providing  a  self  consistent  extrapola¬ 
tion  of  the  data  into  unmeasured  regions.  Since  the  Hilbert 
transform  (1)  is  a  linear  integral  transform  portions  of  ^nr(  uj)  , 
the  transforms  of  which  are  closed-form  analytic  functions,  can 
be  removed  leaving  a  residual  function  whose  transform  must  be 
attacked  numerically.  The  portions  removed  are  chosen  sc  that 
they  smoothly  fit  to  A^(a>)  itself  near  the  boundaries  of  measure¬ 
ment  and  furthermore  give  reasonable  extrapolations  of  ^nr(  m) 
into  the  unmeasured  regions.  Thus  (l)  becomes 

cp(uj)  =  cp0(o>)  +  9(oj)  (5) 


where 


P  r+“  gU1)  dm* 

rr  J  u>’  “  o> 

"  CO 


(10) 


where  g(uj)  is  the  function  which  was  used  to  fit  j&nr(  w)  near 
the  boundaries  of  measurement  and  to  extrapolate  into  the  unknown 
regions.  g(uj)  should  be  chosen  from  among  functions  with  known 
analytic  Hilbert  transforms  «pq(  «))  •  If  vo  d.  note  i\rc{  &)  -  g(  y,) 
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by  A(co)  then 


a(  a1 )  dm1 

to'  “  CB 


(11) 


By  construction  a(u)‘ )  is  a  compact  function  which  is,  in  practice, 
also  bounded.  It  is,  therefore,  square  integrable  and  the 
Hermite  analysis  technique  can  be  used  to  compute  e(®)«  It  is 
then  only  necessary  to  choose  appropriate  g(  «>) .  Our  procedure 
has  been  as  follows.  The  dielectric  constant  or  index  of  refraction 
of  a  system  of  classical  chahge  oscillators  with  given  resonant 
frequencies,  damping  constants,  and  oscillator  strengths  c=a  IS 
well-known.  For  exampleythe  index  of  refraction  can  be  approxi¬ 
mated  (for  a  dilute  gas  of  each  oscillators)  as 

rtfa)  =  1.0  +  E - j— \ -  (12) 

i(«^  -  o>  )  +  ir^® 


in  which  the  are  the  resonant  frequencies,  the  the  damping 
constants  (linewidths)  and  the  essentially  the  oscillator 
"strengths".  The  complex  amplitude  reflection  coefficient  of 
a  bulk  sample  of  material  with  this  ft0is 

P0U)e1<l>0  =  PaU)  =  <fy.)  -  l)/(fij«)  +  1)  (13) 


By  adjusting  the  various  parameters  in  (12)  we  can  find  an  ^(od) 
for  which  the  resulting  Anpo(u))  will  smoothly  and  closely  approximate 
/nr(<j))  near  the  data  boundaries.  The  same  function  (13)  provides 
a  natviVrsui  choice  of  extrapolation  of  r(w)  and  hence  4nr(o>)  into 
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the  unmeasured  regions.  cp0(  u»)  *  the  analytic  portion  we  have 
extracted  from  the  total  Hilbert  transform  is  then  given  by 

*<•>  - tan_1  -57  t  -  <») 

n0^w)  +  0))  -  1 

where,  of  course,  n£ <«)  and  k(< o)  are  the  real  and  imaginary  parts 
of  Aj[uj).  For  the  data  analyzed  in  this  contribution  we  have  found 
that  choosing  two  frequencies  for  classical  oscillators  and 
adjusting  the  other  parameters  r2>  v-^,  v2)  allows  us  to  obtain 

a  sufficiently  smooth  fit  to  the  data  in  the  relevant  regions. 

p  —  h 

Note  that  R(to)  =  1 1*(  a>)  |  has  a  dominant  m  dependence  at  suf¬ 
ficiently  high  frequencies  in  the  extrapolation  here  employed. 

The  entire  data  analysis  problem  has  now  been  reduced  to  the 
computation  of  Hermite  expansion  coefficients  of  a(  o>)  and  matrix 
multiplication  of  these  to  generate  the  expansion  coefficients 
of  e(uj).  The  next  section  is  devoted  to  the  details  of  the  pro¬ 
cedure  we  have  described  as  it  has  been  applied  to  data  on 


SbgSe^  and  BaF2  single  crystals 
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III.  APPLICATION  TO  REAL  SOLIDS 

A.  Optical  Functions  of  Sb2Se3. 

Sb2Se3*  in  crystalline  form,  has  an  orthorhombic  lattice  with  a 
rather  complex  unit  cell  containing  four  units  of  Sb2Se3  with  two  non¬ 
equivalent  sites  both  for  the  antimony  and  selenium  atoms.  The  space 
group  is  (Pnma).  The  optical  functions  are  anisotropic,  however, 
in  this  work  we  ignore  that  fact  as  it  is  not  pertinent  to  our  purpose 
of  illustrating  the  use  of  our  method  of  Kramers-Kronig  analysis.  In 
Fig.  1  is  shown  the  reflectivity  to  unpolarized  light,  in  the  frequency 
range  1.5eV  to  ll.OeV,  of  a  cleavage  plane  of  Sb2Se3  as  measured  by 
Wood6  et  al .  This  data  is  in  close  agreement  with  the  data  of  Shutov® 
et  al  in  overlapping  measured  regions  of  the  spectrum.  The  data,  as  may 
be  seen,  indicates  an  absorption  edge  near  l.OeV  and  a  series  of  peaks 
through  the  high  reflectivity  region  from  the  absorption  edge  through 
the  visible  and  extending  into  the  vacuum  ultra-violet  region.  This  data, 
of  course,  must  be  extrapolated  into  the  unmeasured  regions  by  the 
methods  we  have  outlined  in  Section  II,  above.  The  remainder  of  this 
subsection  details  the  extrapolation  and  other  aspects  of  the  method. 

The  data  were  analyzed  as  follows.  It  was  found  that  a  minimum  of 
two  sets  of  classical  charge  oscillators  with  different  resonant  frequenci 
was  required  in  order  to  extrapolate  the  data  consistently  both  at  high 
frequencies  and  at  low  frequencies  using  (12)  and  (13).  The  criteria  for 
good  extrapolation  that  were  employed  were  (a)  the  residual  portion,  i.e., 
a(u)  in  (11)  should  be  as  small  as  reasonably  possible  at  the  boundaries 
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of  the  data,  and  (b)  a(w)  should  be  rather  smooth  function.  These 
criteria  insure  the  applicability  of  the  analysis  in  II  above  as  well  as 
good  convergence  of  the  Hermite  expansion  of  a(w).  For  the  reflectivity 
spectrum  in  Fig.  1  it  was  found  that  a  good  extrapolation  was  obtaining 
by  choosing  n0(«)  from  (12)  as 


n0  =  !  + 


20 


90 


((2.5)2-u2)  +  i(l)U)  +  (  (11.8)2-u>2)  +  i  (8.5)(u) 


(15) 


Of  course  the  function  above  is  unit-free,  however,  the  units  of  frequency 
in  which  it  can  be  interpreted  are  electron  volts.  Hence  we  have  used 
charged  oscillator  systems  with  resonant  frequencies  2.5  and  11.8  eV. 

Fig.  2  shows  the  function  a(u»)  which  is  obtaining  by  using  for  the  g(u) 
of  section  II  In  p0U)  as  computed  from  (15)  and  subtracting  that 

function  from  In  r  (w)  obtained  from  the  data.  The  contribution  to  the 
phase  from  In  p0(u)  is  obtained  from  (14).  a(u)  has  been  expanded  in  a 
Hermite  series  to  order  250  for  which  the  expansion  coefficients  have  been 
computed  by  Simpson's  rule  with  the  intrgrand  evaluated  at  each  0.05  eV 
interval.  Since  the  data  extends  to  only  12.0  eV  and  the  250th  Hermite 
function  becomes  essentially  zero  beyond  the  "classical  turning  point", 
“top  =^501,  it  is  useful  to  scale  in  frequency  in  order  to  optimize  the 
use  of  the  higher  order  Hermite  functions  in  the  expansion.  Thus  the 
Hermite  expansion  coefficients  of  a(«)  are  given  by 


(16) 


where  Ai(aw)  =  a(w). 
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If  a  =  -v.  1.77  we  have  found  that  the  mean  square  error  between  A (u)  and  its 

Hermite  expansion  is  made  minimal.  The  error  function  x(<*>)  is  plotted  in 

Fig.  3^  and  in  that  figure  <x  (ui)>250  is  given  as  4.54  x  10  .  Thus  the 

-3 

expansion  of  a(u)  is  accurate  to  an  r.m.s.  value  of  order  10  .  Since  we 
know  the  matrix  to  be  compact  the  convergent  expansion  of  +(u)  is 
easily  computed  from  (6).  The  logarithm  of  the  reflectivity  amplitude 
and  its  phase,  +(<■>),  are  shown  in  Fig.  4  for  Sb2Se3* 

The  Kramers-Kronig  analysis  is  now  complete  and  the  optical  functions 
of  the  system  can  easily  be  computed  from  (2)  and  related  equations. 

Fig.  5  shows  the  real  and  imaginary  parts  of  the  index  of  refraction  which 
we  have  obtained  and  Fig.  6  the  real  and  imaginary  parts  of  the  complex 
dielectric  function.  It  may  be  seen  that  the  frequency  positions  of 
critical  points  are  preserved  through  the  transform  process  by  comparing 
Figs.  (5)  and  (6)  with  Fig.  (1). 
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in  the  Mossbauer  spectrum  of  crystalline  Sb2Se^  which 
could  be  fitted  by  postulating  two  equally  populated 
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Sn-Se-As  system  has  been  attributed  by  Shkol’nikov*  to  a 
transition  from  SnSe2  to  SnSe  with  increasing  tin  concen¬ 
tration.  However,  Borisova  et  alz  infer  that  SnSe2  is 
found  only  in  the  glassy  state  and  SnSe  only  in  the  crys¬ 
talline  state.  In  the  present  investigation  of  glasses 
of  Snx(As2Se3)1_x  (with  x  -  2.5,5.0,8.0,10.0  atomic  per¬ 
cent)  were  prepared  by  quenching  the  melts  into  a  brine 
bath  at  200°K.  X-ray  patterns  show  distributions  char¬ 
acteristic  of  a  glass;  further,  the  patterns  are  qualita¬ 
tively  similar  for  all  compositions.  The  Mossbauer  study 
confirms  the  SnSe2-SnSe  transition  with  increasing  tin 
concentration.  The  line  width  and  concentration  depen¬ 
dence  of  the  isomer  shift  indicate  a  high  degree  of 
local  order.  The  question  of  the  presence  of  micro¬ 
crystallinity  will  be  discussed. 


*Work  supported  by  A.R.P.A.  through  A.R.O.,  Durham. 
*E.  V.  Shkol'nikov,  Solid  State  Chemistry  ed.  Z.  N. 
Borisova  (Consultants  Bureau,  N.Y.  1966)  p.  142. 

ZZ.  U.  Borisova  et  al,  Sov.  Phys.  -  Semiconductors  4_, 
443  (1970). 


Submitted  by 

lQ. 

Signature  of  APS  Member 

Clyde  W.  Kimball _ 

Same  name  typewritten 

Physics  Department NIU 
0  Address 

DeKalb,  II 1 iuo i 3  60115 _ _ 


APPENDIX  B 

Abstract  Submitted 


For  the  1971  March  Meeting  of 
The  American  Physical  Society 


29  March  -  1  April  1971 
Date 


I 

* 


Physical  Review 
Analytic  Subject  Index 
Number  47.4 


Bulletin  Subject  Heading  in 
which  Paper  should  be  placed 
Semiconductors  -  Optical 
Properties _ 


Optical  Properties  of  SbjSeo  Single  Crystals. * 
B.  VAN  PELT  and  C.  WOOD,  Northern  Illinois  Univ.--The 
optical  absorption  and  transmission  of  single  crystal 
Sb2Se3  has  been  measured  in  the  region  of  tne  absorp¬ 
tion  edge  with  respect  to  the  principal  crystallogra¬ 
phic  axes.  The  optical  transitions  responsible  for 
the  edge  will  be  discussed. 

*Supported  by  A.R.P.A.  through  Army  Research  Office, 
Durham. 


Submitted  by 


Signature  of  APS  member 


Charles  Wood _ 

Same  name  typewritten 

Department  of  Physics. NIU 
Address 


DeKalb.  Illinois  60115 


Please  place  this  paper  together  with  and  preceeding  the 
accompanying  paper  Photoconductivity  in  Single  Crystals 
of  SbgSej  by  Z.  Hurych,  et  al. 


43 


APPENDIX  B 

Abstract  Submitted 


For  the  1971  March  Meeting  of 
The  American  Physical  Society 


29  March  -  1  April  1971 
Date 


Physical  Review  Bulletin  Subject  Heading  in 

Analytic  Subject  Index  which  Paper  should  be  placed 

Number  47.5  Semiconductors  -  Transport 

Properties _ 


Photoconductivity  in  Single  Crystals  of  SboSe^ . * 
Z.  HURYCH,  C.C.  WANG  and  C.  WOOD,  Northern  Illinois 
Univ. — The  activation  energy  AEa,  the  response  time  t 
and  the  spectral  distribution  of  photocurrent  were 
measured  for  p-type  semiconducting  orthorhombic  single 
crystals  of  Sb2Se^  as  a  function  of  the  orientation  of 
the  electric  vector  EQ  of  the  polarized  light  and  of  the 
drift  field  Ej  with  respect  to  the  three  crystallogra¬ 
phic  axes,  a,  b,  and  c.  Values  of  AEa  were  1.09,  1.10 
and  1.13  eV  for  E0||a,  b,  and  c,  respectively,  in  good 
agreement  with  corresponding  optical  activation  ener¬ 
gies  and  with  published  data  on  photoconductivity  using 
nonpolarized  light.1  Values  of  t  were  in  the  region  300- 
400p  sec  for  hv  =  l.leV,  and  indicate  presence  of  trap¬ 
ping  levels  inside  the  forbidden  gap.  The  spectral 
distribution  of  t  and  as  a  function  of  ^  and  Ed 
orientations  will  be  discussed. 
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Sb  Se  .»  J.C.  SHAFFER,  R.  AFSIIAR,  C.  WOOD,  and  R.  MUEL¬ 
LER,  Northern  Illinois  Univ. — The  optical  functions  of 
thin  films  are  often  difficult  to  obtain  due  to  compli¬ 
cations  arising  from  the  occurrence  of  interference  in 
their  reflection  and  transmission  spectra.  We  present 
and  discuss  ar.  experimental  and  computational  method  for 
circumventing  this  problem.  We  obtain  from  reflectance 
measurements  via  Kramers-Kronig  analysis  the  optical 
functions  of  films,  sufficiently  thick  (>5p)  that  inter¬ 
ference  effects  are  largely  suppressed.  These  optical 
functions  are  employed  to  compute  exactly  the  reflect¬ 
ance  and  transmittance  for  thinner  (M.,00Qft)  film.,  with 
identical  optical  functions.  The  results  are  compared 
with  experimental  reflectance  and  transmission  of 
thinner  films  of  the  same  composition.  These  films  are 
prepared  simultaneously  with  the  thick  films  by  co-evap¬ 
oration.  The  method  is  applied  to  amorphous  films  of 
SbxSey  and  comparison  made  with  the  optical  functions  of 
crystalline  Sb2Se3« 

*Supported  by  A.R.P.A.  through  Army  Researh  Office, 
Durham. 


Submitted  by 


Same  name  typewritten 


Physics  Department,  NTU 
Address 

DeKalb,  Illinois  60115 


52 


APPENDIX  C 


Abstract  Submitted 
For  the  1971  April  Meeting  of 
The  American  Physical  Society 


26  April  -  29  April  1971 
Date 

Physical  Review  Bulletin  Subject  Heading  in 
Analytic  Subject  Index  which  Paper  should  be  placed 
Number  48.3  Amorphous  Solids  _ _ 


Glass-Crystalline  Transitions  in  Amorphous 

Sb.,Se 
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The  g 

selenides  shows  a  strong  dependence  on  composition.  At 
the  composition  Sb2Se3  the  transition  occurs  at  ^180-200°C 
in  agreement  with  the  nuc!  eation  studies  of  Bolotov  et  al. 
The  transition  temperatures  were  determined  by  differen¬ 
tial  thermal  analysis  and  resistivity  versus  temperature 
measurements  and  are  compared. 
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Thermally  Modulated  Reflectivity  of  Sb^Te- . * 

C.  HOOD,  Z.  HURYCH,  and  B.  VAN  PELT,  Northern  Illinois  v 
Univ. — The  thermally  modulated  reflectivity  has  been 
determined  in  the  visible  and  near  IR  region  of  the 
spectrum  for  single  crystals  of  Sb2Te3  with  the  electric 
vector  of  the  radiation  in  the  cleavage  plane  (e|c).  The 
data  show  the  same  critical  points  as  observed  in  the 
optical  transmission  spectrum  of  very  thin  cleaved  speci¬ 
mens  using  a  photon  counter.  These  results  are  compared 
with  the  low  temperature  reflectivity  data  of  Sobolev  - 
et  al.1 
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MOSSBAUER  STUDIES  ON  AMORPHOUS  AND  CRYSTALLINE  ANTIMONY  SELENIDES 

S.  L.  Ruby 
Physics  Division 
Argonne,  Illinois 

and 

L.  R.  Gilbert*  and  C.  Wood* 

Physics  Department 
Northern  Illinois  University 
DeKalb,  Illinois  60115 

Mossbauer  resonance  spectra  at  4°K  have  been  determined  for  a 

range  of  compositions  of  amorphous  SbSex  films,  where  0.5<x<9,  and  were 

compared  with  the  spectrum  of  crystalline  Sb2Se3.  The  crystal  structure 

of  Sb2Se3  is  orthorhombic  Pbnm  (D^Dg5)  consisting  of  infinite  chains 

parallel  to  the  c-axis.  There  are  two  different  equally  populated  Sb 

sites  and  three  different  Se  sites  within  the  crystal  lattice.  The 

strongest  bonds  are  within  the  chains  and  here  the  Sb-Se  bond  distances 

range  from  2.576A  to  2.7778,  while  the  'non-bonded'  Sb-Se  separations 
o  o 

range  from  2.98A  to  3.74A1. 

The  Mossbauer  spectrum  of  a  powdered  sample  of  crystalline  Sb£Se3 
is  shown  twice  in  Fig.  1.  The  solid  line  in  the  upper  curve  is  the 
(unsuccessful)  result  of  a  computer  fit  assuming  a  single  Sb  lattice 
site.  In  the  lower  fit,  the  computer  was  allowed  two  different,  but 

p 

equally  populated  Sb  sites.  Now  x  shrinks  to  1.2  versus  7.3  in  the 

*  This  research  was  supported  by  the  Advanced  Research  Projects  Agency 
of  the  Department  of  Defense  and  was  monitored  by  the  Army  Research 
Office,  Durham,  under  Contract  No.  DA-AR0-D-31 -127-70-G77. 
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summarize,  it  seems  that  for  a  best  fit  a  range  in  e2qQ  and  Sj  are  re¬ 
quired  averaging  near  the  values  given  above. 

The  above  results  show  that  the  short  range  order  in  amorphous 
SbSe  differs  from  that  of  the  crystalline  solid,  and  indicates  a  close 

A 

conformity  in  the  environment,  i.e.,  a  narrow  range  of  coordination 
number  and  Sb-Se  bond  lengths,  of  each  Sb  atom  in  the  amorphous  solid. 
Electron  diffraction  studies  of  Tatarinova2  lend  support  to  this  view¬ 
point.  Intensity  and  radial  distribution  curves  indicate  that  both 
the  coordination  number  and  nearest  neighbor  distances  decrease  in  the 
transition  from  the  crystalline  to  the  amorphous  phase  of  "SbgSe^". 

(The  composition  may  have  been  Se-rich  from  the  method  of  preparation 
described  in  the  article).  In  contrast,  ks^Se^  and  GaAs  were  shown  to 
retain  in  the  amorphous  phase  the  configuration  inherent  in  the  crystal¬ 
line  phase. 

Tne  above  results  resemble  those  from  previous  Mossbauer  studies 
on  glassy  materials  in  which  Fe203  was  dissolved  in  both  silicate  and 
phosphate  glasses''.  The  microenvironments  of  the  Fe  ion  in  each  case 
were  found  to  be  similar  although  the  coordination  was  tetrahedral  in 
the  silicate  and  octahedral  in  the  phosphate  glass.  Despite  the  glassy 
nature  of  the  material  in  these  cases,  as  well  as  in  the  antimony  work 
done  here,  there  does  not  seem  to  be  a  wide  variety  of  microenvironments. 
However,  in  the  amorphous  Sb-Sex  films,  we  have  a  different  microenviron¬ 
ment  than  in  the  crystalline  material. 

These  results  are  of  significance  to  questions  of  whether  a  vapor- 
quenched  solid  is  truly  a  glass  or  consists  of  micro-crystals;  particularly 
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upper  fit.  The  Mossbauer  isomer  shifts  for  the  two  peaks  were  12.6 
tnm/sec  and  15.5  mm/sec. 

The  amorphous  thick  (^  30  microns)  films  were  prepared  by  vacuum 
evaporation  onto  room-temperature  quartz  substrates  from  large  sources 
of  various  ratios  of  Sb  and  Se.  The  films  were  determined  to  be  amor¬ 
phous  by  x-ray  diffraction.  The  composition  as  deposited  (invariably 
different  from  the  source  material)  was  determined  by  electron  micro- 
probe  and/or  x-ray  emission  spectroscopy.  For  example,  a  source  at  a 
temperature  of  720°C  with  x  *  1.5  (i.e.,  Sb^Se^)  produced  a  film  with 
xr*9:  Amorphous  films  of  various  compositions  in  the  range  0.5<x«9 
were  preapred  and  measured  and  the  results  are  shown  in  Table  1.  In  all 
compositions  of  amorphous  SbSex  only  a  single  peak  was  obtained  which 
was  located  at  an  intermediate  position  between  the  double  peaks  in  the 

crystalline  material,  as  Illustrated  in  Fig.  2. 

121 

In  Sb  Mossbauer  spectroscopy,  working  between  nuclear  levels 
with  spins  of  7/2  and  5/2,  a  pure  quadrupole  Interaction  gives  rise 
to  an  asymmetric  octet.  In  the  present  case,  the  quadrupole  Inter¬ 
action  is  not  large  enough  compared  to  the  minimum  line  width  of  7.1 
mm/sec  to  resolve  these  individual  lines.  It  is  found  (see  Fig.  1) 
that  the  data  for  the  amorphous  films  can  be  moderately  well  fitted 
by  a  single  site  where  each  member  of  the  octet  has  a  normal  width 
r=2.8  mm/sec,  the  value  of  e^qQ  Is  10.5  ±  1  mm/sec  and  the  isomer 
shift  is  -13.5  ±  0.5  mm/sec.  By  adding  more  freedom  to  the  computer 
program  (allowing  more  sites),  can  be  further  decreased.  To 
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in  systems  where  the  glass  cannot  be  readily  formed  from  a  melt-quench, 
as  in  the  case  of  Sb^e^.  We  can  state  that  if  the  film  is  an  assemblage 
of  microcrystals,  the  crystal  structure  does  not  correspond  to  the 
equilibrium  phase  and  the  more  probable  conclusion  is  that  it  is  amor¬ 
phous.  The  occurrence  of  a  softening  temperature  (T^)  in  the  films 
would  support  this  conclusion. 

The  authors  wish  to  thank  Dr.  Kimball  for  many  helpful  discussions. 
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TABLE  1 


SbSe 

X 

Phase 

X2 

SI 

e2qQ 

e 

msm 

Amorphous 

.374 

-13.97 

9.64 

3.39 

■9 

Amorphous 

.575 

-14.42 

9.70 

2.68 

2.0 

Amorphous 

.531 

-113.38 

10.85 

3.26 

1.5 

Amorphous 

1.206 

-13.65 

10.90 

3.04 

1.5 

Cryst. 

.253 

-12.60 

7.86 

2.48 

-15.50 

8.40 

2.49 

0.5 

Amorphous 

.466 

-13.00 

11.50 

2.58 

BO 
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The  Optical  Properties  of  Sh2Se3-Sb1Te31) 

By 

W.  Procarione2)  and  C.  Wood 

The  optical  properties  of  single  crystal  specimens  of  Sb2Tc3-SbjRej  solid  solutions  have 
been  investigated  in  the  photon  energy  range  0.04  cV  to  0.5  cV.  A  wide  range  of  solid 
solutions  having  hexagonal  structure  were  found  to  exhibit  pro|>ertics  closely  corresi>onding 
to  the  isostructural  compounds  Ui2Tc3 and  Bi2Scs.  Small  indirect  band  gaps  approaching  se- 
mimetallic  values  were  obtained.  The  high  free- carrier  concentrations  gave  rise  to  plasma 
effects  in  the  near  infra-red.  In  contrast  Sb2Se2  is  orthorhombic  with  a  direct  band  gap 
near  1  eV.  The  chemical  bonding  concepts  applied  to  similar  systems  seem  applicable  with 
only  slight  modification. 

Dio  optischcn  Eigcnsehaften  von  Sb2Te2  -Sb2Nc2-Misohkristallen  wurden  im  Bereich  von 
0,04  eV  bis  6,5eV  untcrsucht.  Die  Festkorperlosungen  haben  innr.rhalb  weiter  Grenzen  die 
gleiche  hexngonale  Struktur  wic  die  isostruktnrcllrn  Verhindungcn  Bi2Te3  und  Bi2Se2  und 
zeigen  sehr  verwundte  optischc  Eigenschaftcn.  Kleine  indiiukte  Encrgieschwellen  ahn- 
lich  denen  in  Halblcitcrn  wurden  bcobachtct.  Die  groBc  Konzentralion  von  f'cicn  Ladungs- 
tragern  verursachte  Plasniaersclieinnngcn  im  nahen  Ultrarot.  Sb2Se3  ist  orthorho.nbisch 
mit  einent  direkten  verbotenen  Encrgicband  von  1  eV.  Die  chemischen  Bindungskonzcpte 
konnen  walirscheinlicli  mit  uur  geriugen  Veranderungcn  auf  ahnlichc  Systcme  angewendet 
werden. 

1.  Introduction 

Binary  compounds  of  the  form  A'B'1  have  liecn  studietl  extensively  in  the 
past.  Solid  solutions  of  these  compounds,  particularly  the  isostructural  systems 
of  Bi.Tej-BijjSc;,  and  Bi2Te3-.Sb2Te3,  have  received  considerable  attention  in 
the  literature.  A  few  serious  attem|>ts  have  been  made  to  study  the  detailed 
band  structure  of  these  compounds  and  their  solid  solutions;  noteworthy  among 
these  being  the  optical  work  of  Greenaway  and  Harbeke  [lj  on  BijTey-BijSe,, 
and  Sobolev  et  al.  (2.  3],  on  Hi2Te3,  Bi,Se3.  and  Sb2Tc3. 

Except  for  some  thermoelectric  datu  by  Teramoto  and  Takayu..agi  [4J,  the 
information  on  the  physical  properties  of  solid  solutions  of  Sb2Te.,  and  Sb^e^ 
appears  to  be  almost  non-existant ,  However,  solid  solutions  of  these  two 
compounds  represent  an  interesting  system  to  study.  Although  possessing 
different  ciy8tallographie  structures,  Sb2Se3  is  orthorhombic  (DJJ  D»)  and 
Sb2Te3  is  rhonihohedral  (CjjDJ.i)  or  hexagonal  on  a  larger  unit  cell,  a  consider¬ 
able  range  of  solid  solubilities  exist  in  (his  system.  X-ray  analysis  byTeramoto 
and  Takayanagi  1 4]  revealed  (hat  a  single-phase  hexagonal  structure  exists 
between  Sb»Te3  and  Sb/l'cSe.  with  the  hit  tic;  constant,  a,  decreasing  approxi¬ 
mately  from  4.27  A  to 4. 12  A  and  the  lattice  constant,  c,  deereasing  approximately 
from  .'10.4  A  to  251. ft  A  with  increasing  Sb2Sc3  content.  Higher  concentrations 
of  SIdScj  resulted  in  an  admixture  of  hexagonal  and  orthorhombic  phase. 

')  Partially  supported  by  A.B.P.A.  through  the  Army  Hesenrch  Office,  Durham, 

•)  This  work  partially  fulfilled  the  requirements  for  the  M  S.  Degree  in  Physics. 
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As  a  start  toward  understanding  the  physical  properties  of  this  system,  and 
obtaining  a  knowledge  of  the  band  structure,  rvfleetior  and  transmission  meas¬ 
urements  on  single  crystals  in  the  range  0.04  eV  to  6.5  eV  were  investigated. 
Because  of  the  marked  directional  preference  in  eleavage,  (0001)  for  hexagonal 
and  (010)  for  orthorhombic  crystals,  and  the  difficulty  of  preparing  good 
reflecting  surfaces  perpendicular  to  the  cleavage  planes,  optical  measurements 
have  been  made  only  with  the  electric  vector  in  the  eleavage  plane. 

Comparison  o!  this  data  with  the  results  of  Greenaway  and  Harbeke  [1]  on 
Bi,Tc3-BiJSe1,  and  Shutov  et  al.  [3]  on  t  he  end  components,  SbjTe,  and  Sb^Se, 
has  shown  many  similarities  between  these  systems  but  has  raised  questions 
concerning  the  chemical  bonding  concepts  formulated  to  explain  the  properties 
of  this  class  of  materials. 


2.  Experimental  Procedures 
2.1  PrrpmrmlioH  of  nmterMo 

All  starting  materials  were  of  six  D’s  purity.  The  antimony  and  tellurium 
were  obtained  from  COM1NCO,  the  selenium  from  Norddcutsche  Raffinerie. 
Further  purification  was  carried  out  on  those  elements  used  for  preparing  the 
end  components  of  the  system,  Sb/’V,  and  Sl^So,.  Fni  these,  the  antimony  and 
tellurium  were,  nut  gassed  and  sealed -oil  under  a  vacuum,  and  then  zone-refined 
in  a  carbon  crucible.  The  selenium  was  out  gassed  and  sublimed  under  vacuo. 

Single  crystals  or  large  grain  size  ingots  were  grown  by  the  Bridgman  tech¬ 
nique  or  by  a  horizont.il  zone-refiner  with  background  beating,  described  else¬ 
where  (3).  Specimens  for  optical  mer.  .nmnents  were  cleaved  from  these  ingots 
down  to  a  thickness  of  a  few  tens  of  microns  hy  (lealing  a  layer  from  the  surface 
using  Scotch-tapc,  or  by  mounting  a  thick  s|iecimcn  with  epoxy  onto  an  alu¬ 
minum  plate  |K>ssessing  a  small  aperture  and  then  pulling  the  bulk  of  the 
material  off  the  plate,  leaving  a  thin  layer  lichind. 

Chemical  analyses  were  perfn rnicd  nil  sections  of  these  ingots  to  test  for 
uniformity  of  composition  and  approximate  corrcs|ioiidcr.rc  with  the  nominal 
composition  was  obtained. 

22i  Op lirml  mrasHrrinmts 

Transmission  and  reflectivity  measurements  in  the  energy  range  6.5  eV  to 
2.0  eV  ^wavelength  0.18  pm  to  2.6  pm)  were  performed  in  a  Cary  double-beam 
14R  spectrophotometer.  A  Cary  attachment,  based  on  n  V-W  configuration  by 
Strong  |G)  was  used  to  measure  reflectivity  at  near  normal  (<  12°)  incidence, 
using  an  aluminum  mirror  as  a  standard  reference.  This  attachment  requires 
two  reflecting  samples. 

A  Beckman  fit  12  spectrophotometer  was  used  for  transmission  and  reflectiv¬ 
ity  measurements  in  the  range  2.0  eV  to  0.04  eV  (2.3  pm  to  33  pm).  A  standard 
Beckman  attachment  was  used  to  measure  reflectivity  at  an  angle  of  incidence 
of  ss2S°.  A  long  wavelength  limit  of  35  pm  was  imposed  on  the  reflectivity 
measurements  hy  the  CsISr.  lenses  in  tlu  iillaeliuicnt. 

3.  Kefleetivily 

The  room-tempcratiire  reflectances  (AM  of  the  SI i.T'-.,  S1);S»‘3  solid  solntions 
in  the  range  0.04  eV  to  2.0  eV  are  shown  in  Fig.  I.  The  pronounced  minima 
appear  to  be  associated  with  the  charge  carrier  plasma  frequency;  in  this  ease, 
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FI*.  1.  The  long- wavelength  reflectance 
Spectra  of  8h1Tet-Sli1Sei  solid  anluliona. 
The  numlxra  refer  to  at°„  Sli,Tf, 


(nerg/feV) 


a  hole  plasma,  since  all  materials  were  determined  to  be  p-type  by  sign  of  the 
Seebeck  coefficient.  A  calculation  of  the  approximate  position  of  wp,  the  plasma 
frequency  minimum  [7]  for  Sb2Te3  from 


p  r‘ 
in*  f„  ( n  'i 


a  e 

l)~  m*fie0(n'i  -  I; 


(1) 


where  p  charge  earlier  concentration,  r  electronic  charge,  in*  effective  mass, 
e0  permittivity  of  free  space,  a  conductivity,  /<  mobility,  and  wr  the  refractive 
index  in  the  non-dispersive  region  of  the  spectrum,  using  values3)  of  p,  m*,  and 
»c  from  Testardi  et.  al.  [8|  gave  a  value  of  7  ion  w  hich  is  in  reasonable  agreement 
with  experiment.  Furthermore,  the  plasma  frequency  minima  for  ShgTe,  samples, 
doped  w  ith  various  amounts  of  copper,  shift  w  ith  the  square  root  of  conductivity 
giving  additional  confirmation  of  this  assignment. 

The  minima  for  the  SbJTe3-Sb,Sc3  solid  solutions  move  toward  longer  wave¬ 
lengths  with  increasing  amounts  of  SloSo.,.  Again,  the  shift  in  position  corre¬ 
lates  surprisingly  well  with  the  square  root  of  the  conductivity  considering  no 
allowance  was  made  for  mobility  variations. 

In  Fig.  2  the  reflectance  values  at  higher  energies  are  shown  and  pronounced 
maxima  are  observed  in  tic  range  I  cV  to  (5  eV.  By  analogy  with  BigTe3,  these 
maxima  are  probably  due  to  direct  iuteiband  transitions  which  should  be  as¬ 
sociated  with  critical  points  at  or  between  the  /V000;>  and  Z(\  1 1>  points  of  the 
Brillonin  zone  1 1  |.  It  this  band  model  is  correct,  then  aiw  line  structure  in  the 
curves  would  be  due  to  a  spin  orbit  split  valence  band.  However,  this  model 


:l)  Published  values  by  vai  ions  investigators  of  p  (mill  a)  for  iimloped  Sb2Te3  are  in  close 
agreeiiient.  The  approximate  value  of  w,.  is  oluaiueil  from  the  curve  of  n -  —  k*  versus  A*  in. 
the  free  earner  absorption  range  by  extrapolating  back  to  zero  wavelength. 
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-  fnergyleV) 


Fig.  3.  The  reflectance  and  trarumlH»lon  spectrum  of 
8UKc, 


Fig.  2.  The  «hort- wavelength  reflectance  spectra  of 
8l»*Te,-8l>,8c,  soli<1  aoluUons 


has  been  rritieiscu  by  Brust  [9]  on  the  grounds  that  the  assignment  of  all  peaks 
to  the  higli  symmetry  p Z  points  is  too  restrietive.  Furthermore,  Sobolev  et  al. 
[2]  found  no  evidence  of  a  triplet  structure  in  their  reflectivity  data  on  Bi/Te*, 
BijScj  and  Sb/Tcj  and  concluded  that  the  reflectivity  maxima  have  no  spin-orbit 
nature.  The  band  model  of  Borghesc  and  Donato  [10]  was  also  too  restricted 
to  high  symmetry  points  to  satisfactorily  explain  their  data  and  they  suggested 
that  assignments  to  the  A(D)  <100)  points  should  not  be  cxludcd  if  the  triplet 
interpretation  is  omitted.  It  apjrears  further  calculations  are  needed  to  unravel 
the  band  structure  of  this  class  ot  materials. 

The  reflectance  data  for  Sb,Sos  is  shown  in  Fig.  3.  The  pronounced  minimum 
in  the  curve  is  associated  with  the  closest  approach  of  the  valence  and  conduction 
bands  and  corresponds  to  a  direct  transition. 

4.  Absorption 

Transmission  measurements  were  obtained  only  on  those  compositions 
containing  the  higher  concentrations  of  Sb^Sca  and  are  shown  in  Fig.  4.  For  the 
high  Sb,Tej  concentrations,  no  measurcable  transmission  was  obtained  even 
for  samples  less  than  10  pm  thick. 

The  absorption  coefficient.  («)  and  hence  the  extinction  coefficient  (k  =  a  A/4  ji) 
was  calculated  from  the  ratio  of  transmitted  (I)  to  incident  intensity  (J0) 

7  (1  -*)*c— 

/„  1  -  7?*e-2«*  ’  V  1 

where  *  is  the  sample  thickness,  assuming  no  multiple  interference  effects. 
Calculations  showed  this  equation  to  be  valid  even  for  those  eases  where  inter- 
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Fig.  4.  Tile  transmission  spectrum  of  Sb.Tej-Sb.Re,  FI*.  5.  The  vsrlstlon  of  the  product  of  absorption  eo- 
solid  solutions  efficient  (a)  and  refractive  index  (n)lntlio  free  carrier 

absorption  region  with  the  square  of  the  wavelength 
in  Sb/Tci-St^Se,  solid  solutions 


ference  fringes  were  observed,  leading  to  errors  in  a  of  less  than  2%  for  these 
samples.  Having  obtained  k ,  the  refractive  index  (n)  was  then  determined  from 


n  = 


1  +  * 
1  -  R 


.  ]PJl1 _ fci 

\  (1  -  R)*  ' 


(3) 


Values  so  obtained  agreed  very  closely  with  the  values  of  n  determined 
from  the  positions  of  the  interference  fringes.  The  product  n  a  in  the  free  carrier 
absorption  range  is  plotted  against  A*  in  Fig.  5  and  is  shown  to  obey  the  relation 

[7] 


na  = 


4  71* 


A*  es  p 

e*  £„  ml  nv ' 


(4) 


The  square  root  of  the  absorption  coefficient  was  plotted  against  energy 
(Fig.  0)  for  each  composition  and  found  to  yield  two  straight  lines  conforming 
to  the  equation  [11] 


a 


A 


~(hv  -  Es-  kO)2  (A  v  -  E,  +  k  Of 

!-«*(-#)  •*(£)-> 


(5) 


Fig.  0,  Square  root  of  the  alworption  coefficient 
vpmiis  incident  photon  energy 


~[nergy(eV) 
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for  the  simplest  ease  of  indirect  allowed  transitions,  where  k  6  is  the  single 
phonon  energy  required  to  conserve  momentum.  By  analogy  with  Bi,Te,  [1] 
these  indiieet  transitions  should  occur  somewhere  along  the  PA  or  P  D  direc¬ 
tions  of  the  Hrillonin  zone.  As  is  generally  the  case  for  if(  type  critical  points  no 
structure  in  the  room-temperature  reflectivity  curves  corresponding  to  these 
indirect  transitions  was  observed  [12]. 

An  arbitrary  selection  of  an  absorption  coefficient  value  for  the  absorption 
edge  [1]  did  not  show  any  regular  trend  in  energy  gap  with  composition.  Values 
of  the  energy  gap  ( A'„)  and  k  0  determined  from  the  intercepts  of  the  curves  with 
the  abscissa  in  Fig.  6  are  listed  in  Table  1.  The  Et- values  show  a  decreasing 
trend  with  increasing  SbjSe,  content. 

The  exact  values  of  these  energy  gaps  are  uncertain  because  of  the  indication 
of  considerable,  charge  carrier  degeneracy  from  Fermi-level  positions  determined 
by  Secbeck  coefficient  measurements.  The  possibility  of  the  SbjTej-SbjStj 
solid  solutions  having  light  and  heavy  hole  valence  bands  as  in  SbxTca  [13,  14], 
or  that  both  electrons  and  holes  are  contributing  to  the  Seebeck  coefficient  in 
these  low  band  gap  materials,  does  not  allow  an  accurate  designation  to  be  made 
of  the  position  of  the  Fermi -level  from  Seebeck  coefficient  measurements. 
Moreover,  corrections  for  degeneracy,  as  carried  out  by  Sehr  and  Testardi  [15] 
or  by  Greenaway  and  Ilarhckc  [1],  did  not  allow  for  the  fact  that  optical  transi¬ 
tions  may  take  place  from  levels  approximately  4  kT  (equivalent,  to  0.1  eV 
at  rooni-tciniM'i aturc)  away  from  the  Fcnni-lovol  [16],  i.e.  at  an  energy  spread 
greater  than  the  correction  applied.  Therefore,  no  reliable  estimate  of  the  Bur- 
stcin  shift  [16]  can  be  made. 

For  pure  SbjScj  the  energy  gap  occurs  at  about  1  eV  and  in  contrast  to  the 
statement  of  Shutov  et  al.  [3],  appears  to  arise  from  a  direct  allowed  transition 
[17],  K  xceilcut  correspondence  is  seen  in  Fig.  3  between  the  reflectivity  mini¬ 
mum  and  the  point  of  maxim  urn  slope  of  the  absorption  edge  [7J. 


Table  1 


Sample 

No. 

Nominal  | 
%  SbjTe,  | 

1 

Resistivity 

(Oem) 

(300  °K) 

Seobeek  eoc.ff. 
(pV/deg) 
(300  °K) 

Fermi 
level  **) 
(cV) 

i 

I  (eV)  | 

1  1 

6 

(°K) 

Amin 

(pm) 

1 

100 

2.18x10* 

!H)  ±  2 

0.07 

0.21  *) 

400  to 

8.5 

1 

1000*) 

o 

no 

4.18x10-* 

45  ±  3 

0.14 

— 

— 

10 

1 

78 

— 

75  ±  3 

0.09 

— 

— 

11 

4 

(17 

4.27  X  I0-* 

43  ±  4 

0.17 

0.10 

775 

13 

5 

00 

4.83x10* 

48  ±  3 

0.14 

0.14 

600 

15 

0 

no 

5.00  x  I0-* 

97  ±  2 

0.00 

0.09 

000 

19 

7 

41) 

8.00  x  fit-* 

01  ±  4 

0.11 

0.04 

520 

23 

8 

33 

11.28x10** 

(50  +  2 

0.10 

0 

500 

24 

*)  Values  obtained  from  f  If* ). 

**)  Fermi  level  values  ealculiitcd  from  llie  Neelicvi;  coefficient  assuming  a  oif-carrier 
single  vai  -nee  Imiul  model. 

5.  Discussion 

In  view  of  the  fact  lhat  Sb2Tc3  is  isost met nml  with  Bi,Tc3  and  BijSej,  the 
highest,  valence  hand  and  lowest  conduction  band  would  be  expected  to  consist 
of  six  ellipsoidal  isoenergetie  surfaces  |2|  and  the  Brillouin  zone  would  be  ex- 
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pected  to  exhibit  the  same  symmetry  as  Ui.Te,  [1],  The  similarity  in  optical 
properties  is  therefore  not  surprising.  SbjSe,  differs  from  Sb,Te,  in  its  crystal 
structure.  Thus,  some  differences  in  optical  properties  would  be  expected  and 
are  found,  e.g.  the  smallest  energy  gap  in  Sh^Se,  is  a  direct  one,  and  splitting  of 
several  SbjSe,  reflectivity  peaks  (A«.j)  occur  due  to  lower  crystal  symmetry  [3]. 

Drabble  and  Goodman  [18]  have  proposed  a  layer  structure  for  the  Bi/Te, 
class  of  compounds,  each  layer  consisting  of  similar  atoms  stacked  in  a  sequence 
TeO)-Bi-Te<2,-Bi-Te(|)  perpendicular  to  the  three-fold  axis.  The  superscripts 
denote  differently  bonded  tellurium  atoms:  the  Tc(2>  atoms  being  Burrounded 
almost  octohedrally  by  covalently  bonded  Bi  atoms,  with  a  weak  ionic  compo¬ 
nent;  the  Tc(l>  atoms  have  three  nearest  neighbor  covalently  bonded  Bi  atoms, 
with  a  strong  ionic  component,  and  three  van  del-  Waal  bonded  Tc<'>  atoms  in 
the  adjacent  layer.  This  model  was  used  by  Drabble  and  Goodman  to  explain 
the  results  of  Austin  and  Slieard  [19]  on  the  variation  of  the  optical  energy 
gap  of  BigTcy-rScx  alloys,  and  was  used  later  by  Greenaway  and  llarbckc  [1] 
to  explain  their  data  on  similar  alloys.  It  was  promised  that  if  tellurium  is 
replaced  by  selenium  then  the  more  electronegative  selenium  atoms  will  sub¬ 
stitute  preferentially  for  the  more  weakly  bound  Te*2>  atoms.  The  increased 
bond  strength  will  eansc  an  increase  in  the  energy  gap.  Eventually,  when  all 
Te(2)  sites  are  occupied,  selenium  then  substitutes  for  Tc(|)  atoms,  attracting 
charge  away  from  the  Bi -tSe(2)  bonds,  thuB  causing  the  energy  gap  to  decrease 
again. 

It'  appears  more  likely  that  selenium  is  substituting  for  the  T'(|>  atoms  in 
the  Sb^Tej-BbjSej  alloys.  The  general  trend  ol  decreasing  band  gaps  and  the 
transition  from  hexagonal  to  mixed  phase  at  the  composition  Sb,Tc,Sc2  lend 
support  to  this  view.  The  high  hole  concentration  always  observed  in  SbsTe3  is 
thought  to  be  due  to  a  wrong-atom  defect  [20]  with  a  stoichiometric  excess  of  Sb 
(providing  acceptor  states)  on  Tc  sites.  Although  it  is  probable,  in  view  of  the 
above  model,  that  Te  vacancies  would  occur  preferentially  on  Te(2)  sites  it 
does  not  necessarily  follow,  as  Tcramoto  and  Takayanagi  [4]  supposed,  that  it  is 
energetically  favorable  for  excess  Sb  to  occupy  the  Te(2)  sites  rather  than  the 
Te**>  sites. 

It  appears  that  more  definitive  X-ray  studies  arc  required  to  resolve  these 
problems. 
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THERMALLY  MODULATED  PLASMA  FREQUENCY  OF  FREE  HOLES 
IN  THE  Sb2Te3^SeJC  COMPOUNDS* 
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( Received  3  November  1170  by  R.H.  Silsbee) 


The  plasma  resonance  frequency  oj.  of  free  holes  in  a  single-crystal 
of  the  narrow-band  semiconductor  Sb2  T e 3.x Se*  was  obtained  using 
thermally-modulated  reflectivity.  This  method  gave  accurate  values 
of  even  in  the  case  of  broadening  due  to  strong  free-carrier  and 
interband  absorption. 


SINGLE  crystal  solid  solutions  of  rhombohedral 
Sb2Te2  and  orthorhombic  Sb2Ses  can  be  prepared 
of  composition  Sb2TeMSe*  for  values  of  x  in  the 
range  0-2.'  All  of  these  compounds  possess  simi¬ 
lar  physical  properties:  rhombohedral  ( CnD\d  ) 
crystal  structure;  p- type  conductivity  with  values 
of  10  *  —10  3  (Ocm)  1 ,  decreasing  with  x;  similar 
reflectivity  spectra  with  two  sharp  reflectivity 
windows,  one  in  the  u.v.  region  at  A  ~  2,800  A 
for  all  samples,  and  the  other  in  the  i.r.  region  at 
A~  9  *5/i  for  Sb2Te3  and  at  longer  wavelengths 
with  increasing  x. 

The  long  wavelength  edges  of  these  i.r.  reflec¬ 
tivity  windows  were  associated  with  the  plasma 
resonance  frequency  ojp  of  free  holes  in  the 
valence  band.  For  Sb2Tes  this  assignment  of  a 
plasma  window  was  confirmed1  from  the  position 
of  the  minimum  reflectivity.  However,  absorptive 
processes  produce  a  broadening  of  the  resonance 
pattern  as  well  as  a  non-zero  reflectivity  at  the 
minimum,  thus  complicating  the  exact  determination 
of  o>p  from  the  reflectance  spectrum.  The  rather 
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of  Defense  and  was  monitored  by  the  Army 
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DA-ARO-D-31-124-70-G77. 
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Fig.  la,  b:  The  reflectivity  R  (dashed  curve)  and 
modulated  reflectivity  AR/R  (solid  curve)  in  the 
vicinity  of  the  plasma  resonance  frequency  for 
various  compositions  of  Sb2TeMSe^. 
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Table  1. 


Sb2Te, 

Sb2Te2Se 

S^Te^Se,., 

Sb2TeSe2 

Ci)p  [cm-'  ] 

from  thermal  modulation 

1,040 

645 

435 

395 

o>2  [cm*'  1 

from  thermal  modulation 

1,125 

665 

505 

470 

7.3 

16.7 

4.0 

3.70 

ec 

from  extrapolation  A  =  0 

50* 

28 

42 

39 

V 

0.410 

0.158 

0.860 

0.925 

carlcm-1  J 

from  77 

426 

102 

375 

370 

Wp.Icm'1  ] 

from  extrapolation  £r  =  0 

1,025* 

530 

490 

373 

*  Data  from  reference  7. 


high  concentration  of  free  carriers  (p  ~  10'*— 10*®cm"*) 
affects  not  only  the  real  part  of  the  dielectric 
constant  £r  but  also  the  imaginary  part  e,  through 
the  free-carrier  absorption.  Moreover,  in  the  case 
of  p-type  materials,  an  additional  contribution  to 
£r  and  to  £,  due  to  interband  transition  between 
light  and  heavy  hole  bands  generally  exists,2'3 
which  may  be  enhanced  in  the  case  of  spin-orbit 
splitting  of  the  valence  bands.  Therefore,  thermo- 
modulated  reflectivity  experiments  were  performed 
on  the  Sb2Tea.2.Sez  compounds  in  order  to  obtain 
accurate  values  of  o)p  and  to  confirm  the  plasma 
oscillation  assignment. 


All  samples  could  be  cleaved  easily  along 
natural  cleavage  planes  (0001)  thus  allowing  tlx; 
preparation  of  thin  single  crystal  specimens, 
20—50/4  thick.  With  this  sample  orientation  the 
electric  vector  '£  of  the  incident  radiation  was 
perpendicular  to  c-axis  throughout  our  measuie- 
ments.  Square-wave  current  pulses  with  peak 
values  of  from  500mA  up  to  2.0  A  at  a  frequency 
of  11  c/s  were  passed  through  the  samples.  A 
modified  Beckmann  i.r.-12  Spectrophotometer  was 
•ised,  such  that  the  modulated  11  cps  signal  from 
tin.  thermopile  detector  was  fed  to  a  Keithley  802 
Serie.'  Lock-in  Amplifier  and  then  recorded. 

The  thermoreflectance  spectra  of  the  samples 
are  shown  in  Fig.  1  together  with  the  reflectance 
spectra.  All  of  the  samples  exhibit  the  typical 
feature;  of  a  thermally  modulated  p.'asma  frequency, 


Fig.  2.  The  real  part  of  the  dielectric  constant, 
e,  =  n2  -  k*  and  the  free-carrier  absorption 
(n  *  a)  versus  the  square  of  the  wavelength  (A2) 
for  the  Sb2Te3  zSex  compounds. 

i.e.  two  sharp  peaks  of  opposite  sign.*  One  of  the 
peaks  (u >p)  corresponds  to  the  thermal  modulation 
of  the  plasma  frequency  o)p  (in  these  narrow  gap  semi 
semiconductors,  cop  is  modulated  both  by  change 
of  carrier  concentration,  p,  as  well  as  by  change 
of  the  effective  mass  mj),  and  the  other  peak 
(a;,)  corresponds  to  the  modulation  of  the  broad¬ 
ening  parameter  77  =  o)T/o:p,  (where  (t)T  is  the 
collision  frequency  which  express  phenomeno¬ 
logically  l  he  scattering  of  free  carriers  by  phonons). 
The  shorter  wavelength  side  of  the  thermomodulated 
spectra  shows  additional  structure  which  may  be 
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doe  to  tausitioai  between  the  light  and  heavy 
hole  bands.*’*  Light  and  heavy  hole  hands  have 
been  postulated  to  explain  thermoelectric*  and 
Hall  effect  date*  in  this  class  ai  compounds  usd 
appear  to  indicate  a  heavy  to  light  hole  concen¬ 
tration  of  at  least  two  orders  of  magnitude.  This 
large  concentration  ratio,  even  allowing  for  the 
lower  effective  mass  of  light  holes,  would  suggest 
that  a  plasma  resonance  frequency  due  to  the  light 
holes  would  occur  somehweie  in  the  far  i.r.  region 
of  the  spectrum  [sec  equation  (2)],  and  that  heavy 
holes  are  responsible  for  the  plasma  resonance 
observed. 


From  the  two  frequencies  eop  and  eot,  the 
non-dispetsive  dielectric  constant  of  the  lattice 
eL  and  the  broadening  parameter  tj  can  be  deter¬ 
mined  using  the  relations4 


and  R(cot) 


4  GJ* 


.  (1) 


These  results  ate  shown  in  Table  1  together  with 
calculated  values  of  &>r  and  of  the  optical  dielectric 
constant  e  c,  the  latter  was  determined  from  the 
extrapolation  of  the  spectral  dependence  of  n2-k 1 
in  the  free  carrier  absorption  range  to  A  -  0. 

(Fig.  2;  n  and  k  were  obtained  from  transmission 
and  reflectivity  measurements.1 ) 


due  to  large  discrepancy  between  the  wavelength 
at  which  minimum  reflectivity  occurs  and  the  wave¬ 
length  were  e,  -  n*  -  k*  *1  [  e&has  very  strong 
dependence  on  (&>,/o>a)].  A  similar  discrepancy 
can  be  found  in  case  of  SbaTea.7  This  discrepancy 
can  he  due  either  to  interband  transitions  he  tween 
light  and  heavy  hole  bands  which  contribute  to 
the  real,  er,  and  imaginary,  eit  parts  of  the 
dielectric  constant,  similar  to  the  case  of  p-typr 
germanium,3  or  to  a  high  collision  frequency  <ua 
such  that  the  condition  (a>p/u>T'f  »  1  is  not  valid. 
For  either  case,  the  expression: 


4np 

(eo  +  i  o)T) 


(2) 


for  the  total  dielectric  constant  does  not  apply,  so 
equation  (1)  cannot  be  used  without  modification. 
Existence  of  either  of  the  above  phenomena  is 
supported  by  an  absorption  background  resulting 
in  high  values  of  the  free  carrier  absorption  back¬ 
ground  resulting  in  high  values  of  the  free  carrier 
absorption  (~  10*  enf '  )  as  extrapolated  to  A  =  0. 
However,  the  additional  structure  in  the  thermo- 
modulated  spectra  as  well  as  the  calculation  of 
air  from  the  Hall  modility  data7  indicate  that  the 
effect  of  interband  transitions  is  more  probable. 
These  transitions  are  therefore  also  responsible 
for  the  discrepancy  between  up  end  &>;. 


Extrapolation  of  the  aame  curves  to  n1  -k1  «=  0 
gave  values  of  to’  in  reasonable  agreement  with 
the  values  of  cop  obtained  from  thermal  modulation 
while  values  of  determined  from  equation  (1) 
were  quite  different  from  ec,  which  is  obviously 


The  advantage  of  therm omodula ting  the  reflec¬ 
tion  spectra  is  clearly  indicated  by  the  fact  that 
despite  very  st  ong  broadening  due  to  Competing 
absorption  processes,  cop  can  be  determined  with 
high  accuracy. 
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Die  Plastu3r.;sonanxfrequenz  cjp  von  freien  LOchern  in  Einkristallen 
dea  Schmalband-Halbleiters  SbaTea_aSes  wurde  mit  Hilfe  von  tbermisch 
modulierten  ReflektivitBtamessungen  bestimmt.  Mit  dieser  Methode 
erzielt  man  sogar  dann  genaue  Ergebnisse,  wenn  Verschmierung 
dutch  starke  Absorption  der  freien  Trfiger  auftritt  odcr  wenn  Inter- 
bandabsorption  stattfindet. 
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